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Abstract 

Nematicity is a well known property of liquid crystals and has been recently discussed in the 
context of strongly interacting electrons. An electronic nematic phase has been seen by many 
experiments in certain strongly correlated materials, in particular, in the pseudogap phase generic 
to many hole-doped cuprate superconductors. Recent measurements in high T c superconductors 
has shown even if the lattice is perfectly rotationally symmetric, the ground state can still have 
strongly nematic local properties. Our study of the two-dimensional Hubbard model provides 
strong support of the recent experimental results on local rotational C4 symmetry breaking. The 
variational cluster approach is used here to show the possibility of an electronic nematic state and 
the proximity of the underlying symmetry-breaking ground state within the Hubbard model. We 
identify this nematic phase in the overdoped region and show that the local nematicity decreases 
with increasing electron filling. Our results also indicate that strong Coulomb interaction may 
drive the nematic phase into a phase similar to the stripe structure. The calculated spin (magnetic) 
correlation function in momentum space shows the effects resulting from real-space nematicity. 
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INTRODUCTION 



Strongly correlated electron systems, such as the high temperature superconductors and 
other oxides, exhibit complex phase diagrams. In spite of intensive research during the past 
decade, there is no consensus on a possible mechanism for the superconductivity and certain 
new phases that are found from time to time. Recent measurements of magnetoresistivity, 
resistivity, IR reflectivity and neutron scattering in ruthenates[l|, iron-based superconduc- 
tors {2I, y[ and cuprates J4-8] show an interesting phase which locally breaks the rotational 
point group symmetry (C4) along x and y axes (these directions for cuprates, defined along 
the Cu-0 bond directions in the Cu-0 plane) but preserves the translational symmetry. The 
symmetry of this phase greatly resembles the well-known nematic phase in isotropic liquid 
crystals, and hence the phase is called an electronic nematic state. 

A two-dimensional (2d) electronic nematic phase is usually referred to as a phase that 
spontaneously breaks the symmetry of the underlying Hamiltonian associated with inter- 



changing the x and y axes (defined in a certain plane) of the system [25|. Unlike conventional 
liquid crystals, in which the nematic state is due to the special shape of the molecules or the 
anisotropic interaction between molecules, the origin of an electronic nematic phase is the 
local electronic correlation. When C4 symmetry is reduced to C2 symmetry, the system is 
expected to have a 2-fold degenerate ground state. The two ground states will be randomly 
distributed on the 2d plane and will form locally different domains similar to magnetic do- 
mains in an Ising model. The effects of differently oriented domains are very likely to cancel 
out globally (i.e., on a larger length scale) and give rise to a lattice that appears uniform 
(i.e., without the above symmetry breaking). 

It is technically difficult to detect the local electronic nematic symmetry breaking with 
disordered domains. Many experiments focus on Y Ba 2 Cu^O^ +x because its orthorhombic 
structure plays a role similar to an external field that can align nematic domains along the 
same direction, so that it is easy to measure the macroscopic anisotropic effect. However the 
orthorhombicity also makes it hard to distinguish whether the anisotropic effect is caused by 
the distorted lattice structure or the local electronic correlation. Recent publications % |5] 
have reported a clear and strong anisotropic resistivity even when the orthorhombicity is 
well suppressed and the temperature dependence of the anisotropic effect is different from 



that of a lattice distortion. This implies that the in-plane transport anisotropy cannot be 
explained solely by distorted lattice structures. This result was later supported by a cellular 
dynamical mean field theory (CDMFT)[9] study of the two-dimensional Hubbard model by 
Okamoto et al. [10]. They showed that a strong nematic resistivity could be induced by 
a very small orthorhombic distortion as long as the interaction is strong enough to yield a 
Mott transition. 



Another way to identify local nematicity is by using scanning tunneling spectroscopy 
(STS) with a rather high resolution. The STS can detect directly the electronic density 
of states within a nematic domain. In a recent publication, it is reported that the state- 
of-the-art STS measurements show "intra-unit-cell" electronic nematicity [8| on a perfectly 
symmetric lattice along the x and y axes i.e., a lattice that is not orthorhombic. Thus, 
this measurement provides the most direct evidence that the local nematic effect, in some 
strongly correlated materials, is purely electronic and originates from repulsive electron- 
electron interactions. 



The experimental results from STS detection depict a fascinating picture of disordered 
local nematic states fluctuating throughout the undistorted cooper oxide plane, but, to 
our knowledge, there is no simple theoretical model supporting this idea. The motivation 
for this report is the possible existence of a ground state of the one-band Hubbard lattice 
Hamiltonian with spontaneous x — y symmetry breaking. In order to identify and address 
such a locally nematic ground state, a real-space quantum cluster method which can 
treat both the local and global effects at the same time is probably the best option for the 
problem. In this report, the variational cluster approximation (VCA) [12] calculations are 
performed on a ^-symmetric, two-dimensional Hubbard square lattice. The results from 
the VCA calculations are used to make further analysis of this phase. 



The report is organized as follows. The models and methods used in this paper, including 
a brief introduction to the VCA, are presented in Sec. „ Sec. n first examines the existence of 
a local nematic phase and its coexistence with antiferromagnetism (Sec. n and Sec. „). A local 
order parameter is introduced in Sec. „ and the results on the dynamical structure factor are 
discussed in Sec. n while concluding remarks are given in Sec. „ 
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METHODOLOGY 



Hamiltonian 



To explore a possible mechanism for local nematicity, we consider here the simplest one- 
band Hubbard Hamiltonian if on a square lattice 

H = (-tCra Cr'a + H.C.) + U ^ fl^K^ i 1 ) 

(r,r')a r 

where c ra is the annihilation operator for electrons at site r with spin projection a and 
U is the on-site screened Coulomb repulsion, (r, r') denotes summation over the nearest 
neighbors. The parameter t is the overlap integral between nearest neighbors. The two- 
dimensional square (t r , r+x = t r ^ r+y = t) lattice Hamiltonian is symmetric along x and 
y directions, i.e., it has the underlying C4 symmetry. All the energies reported here are 
measured in units of t. 

Variational Cluster Approximation (VCA) 



There are many different numerical methods designed to solve the Hubbard model ap- 



proximately. We wi 
quantum cluster 



1 use the variational cluster approximation in this article. This is a 

n 

11| extension of Potthoff's self-energy functional approach (SFA) [13]. 
The SFA develops a variational process by using self-energy of a "reference system" that 
can be exactly solved to obtain the approximate ground state properties of the "real system" 



in the thermodynamic limit. In Ref. 



13(, Potthoff proved that the grand potential Q of a 



system, with one-particle parameter t and on-site interaction parameter U, can be written 
as a functional of the self-energy X provided the system is not at the critical point of a phase 
transition 13] : 

= TrlnUGo 1 - S)- 1 ) + F[E] 
with = G[S]- X = GofS]- 1 - S 



(2) 



where Go is non- interacting Green's function of the original model and i* 1 [£(£)] is the Leg- 
endre transform of the Luttinger-Ward functional [ijj]. The equations indicate that 
is stationary at the exact self-energy and its value is the exact grand potential of the sys- 
tem. A "reference system" is introduced with the same on-site interaction U but different 
one-particle parameter t', and the reference system can be exactly solved. In the VCA, the 
reference system consists of isolated, identical clusters from a decomposition of the original 
lattice. Because of the universality of the Luttinger-Ward functional with the same U, the 
functional dependence of F[S(t')] will be the same for the original system and the reference 
system. Therefore, F[S] can be extracted from the reference system and used in the original 
system, so that F[S\ can be eliminated. 

ftt[E] = f2 t /[S] - Trln(G Q l - S) + Trln^ 1 - £), (3) 

and the stationary point of the equation is at S = S[t]. However, when evaluating the 
functional Ot[S] for the stationary point, one can only vary the functional within a subspace 
of self-energy S(t') instead of the full space This is where the approximation kicks in. 

By evaluating the stationary points of the functional in the restricted self-energy space, the 
grand potential and the Green's function of the original system can be found approximately. 
It is possible that the approximation comes with multiple stationary points; it is also possible 
that some stationary points will shift or even disappear in a different self-energy space 
because of the approximation, but those not disappearing should be physically meaningful 



13|- 



solutions and the true ground state resides at the state with the lowest grand potential 
The VCA 12] is constructed upon the SFA by choosing the reference system an ensemble 



of independent, identical clusters from a decomposition of the original lattice. The cluster 



can be solved by either exact diagonalization 18j or Lanczos method 26[. The exact eigen 



values and eigenstates of the cluster lead to an accurate implementation of the correlation 
effect, at least within the dimensions of the cluster size, and long range correlation acts like 
a mean- field in this approach. Evaluating the stationary point (s) of f2[S] by varying the 
one particle hopping parameter in the cluster t' provides a good approximation to the exact 
ground state of the original system in the thermodynamic limit. The electronic density n in 
the lattice can be controlled by the chemical potential of the lattice \x. Due to the mean-field 
treatment beyond the cluster dimension, number of electrons per cluster can be non-integral. 
Hence, although the number of electrons in the reference system (the finite cluster) can only 
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FIG. 1. (a) The reference system used in this study: the whole lattice is divided into identical 
2x2 (square symmetry) clusters with the same coupling t' between the nearest atoms which is 
considered a variational parameter. Additional variational parameters are chemical potential // , 
effective Weiss field M and fictitious deviation 5 (see Sec. „ for details), (b) The full square is the 
first Brillouin zone (FBZ); shaded square denotes the reduced, first Brillouin zone (RFBZ) of the 
cluster superlattice. In a cluster representation, in the real lattice, upper-case letters represent 
lattice vectors within a cluster and the tilded letters denote a vector connecting the origins of two 
lattices; in the reciprocal lattice, we use upper-case letters for a vector between origins of two 
RFBZs and the tilded letters to represent the momentum vector within a RFBZ. For example, 
x = X + x in (a) and k = K + k in (b). We use the above definitions throughout this article. 

be an integer, it is indeed possible to achieve any electronic density from the VCA approach. 
Reference System for Solving Nematicity 



In this report, we introduce nematicity locally in the smallest cluster cell; 2x2 (square 
symmetry) clusters are used as the reference system (Fig. QJa)) to tile the infinite square 



and t' is the hopping parameter in the reference system. Actua 



m 



ly. 



one can also apply the 



20| | . The clusters use open 



larger-sized Betts clusters to reduce the size and edge effects 

boundary conditions according to Ref. 12| and t' is a variational parameter for the VCA 
calculation. 

Without coupling among the clusters, it is difficult to directly control the electronic 
density since the clusters have an integral number of electrons which is conserved. In the 



VCA, the chemical potential /x is varied to achieve different electronic density. However, the 
chemical potential of the clusters in the reference system, //, should be different from \x in 
the original system, so // is also treated as a variational parameter. This treatment also 
ensures the thermodynamic consistency of the VCA, i.e., the average electronic density (n) 
calculated from the trace of the Green's function, TrG, and the negative derivative of the 



grand potential —dVl/d/j, are supposed to be the same [21]. 

The most significant improvement of the VCA over cluster perturbation theory |22| is 
that the VCA can describe the spontaneous symmetry-breaking orders by introducing a 
fictitious external field. In this paper, we will discuss two kinds of broken symmetries: 
antiferromagnetic (AF) and nematic. As for the AF order, one more term can be included in 
the reference system to represent the broken AF symmetry: Hm = M e l ® AF (n r ^ — n r ±), 

r 

where Qaf is the AF wavevector and M is the fictitious AF Weiss field whose value is 
determined by the variational procedure 23]. Here we only consider the simplest possibility 
for the AF case, which has a static AF order with Qaf = (7r, 7r). The nematic state gives 
rise to a fictitious difference between hopping along x and y directions; this difference can 
be addressed by introducing an additional term H§ = SY^t'(c^c r+x — c^c r+y ) + H.c, where 

r 

S denotes the fictitious deviation of the reference system from C4 symmetry. The subscript 
r + x(y) indicates the nearest neighboring site along the x(y) direction in the reference 
system. Here, we confine ourselves to 5 > 0, which means that the hopping along the 
x direction of the reference system is larger or equal to that along the y direction. This 
will provide one ground state, and because the ground state for nematic order is two-fold 
degenerate, the other ground state is symmetric to the first ground state with a larger 
hopping along the y direction. Both M and 6 are zero for the symmetric states and non- 
zero for the symmetry-breaking orders. It is important to note that M and 5 are both 
fictitious parameters, which indicate the influence of long-range order on the small clusters, 
so the bigger the size of the cluster, the smaller the fictitious parameters are. This rule can 
be used to justify the validity of the fictitious parameters: the case for M is proved following 
this rule in Ref. {24 1. 

Finally, with these one-particle variational parameters, the reference system is represented 
by a Hamiltonian H' given by 

H>= ]T ( - t'c+ a K<o + H.c.) - n 

(r,r')a r 



+8 Y,[t'(^c r+X - c+c r+y ) + H.c] 

r 

+M e lQAF - R (n rt - n ri ) + n rt n ri (4) 



r 



where H' contains the same interacting term as H in Eq. [T] except for the modified one- 
particle terms which create a trial self-energy for the variational process. A series of trial 
self-energies are determined from the variational calculation which ends at the stationary 
points (which can be maxima, minima or saddle points) of the self-energy functional 
The variational calculation involves four variational parameters (t', /i', M and S) making 
it difficult to find the correct stationary point (s), because it is not possible to search every 
corner of the variational space completely. Therefore, we carry out the search in this report 
as follows: (1) set 5 = and M = to find a trivial solution by the VCA; (2) vary only 
the fictitious deviation 5 and the effective Weiss field parameter M while fixing all the 
other parameters from the previous process; (3) as long as a non-trivial result (5 ^ or 
M ^ 0) shows up, an additional variational process with all the variational parameters 
will be performed to get the final result. The above three-step process is the most efficient 
way we have found in order to obtain a consistent variational behavior within the limited 
computer resources and time. We use a VCA program written by ourselves in FORTRAN. 
The program has been tested extensively and verified with results in some published work 
such as Refs. [l2, 
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RESULTS 



Existence of Local Nematicity 



The first question here is whether a local nematic state really exists in our simple model. 
The general experience from VCA calculations is that as long as the second step of the 
variational process mentioned in Sec. „ gives a nontrivial result, the true stationary point 
obtained in the third step is very likely to be nearby and the nontrivial result usually survives 
in the third step. Therefore, we can use results from the second step to test the validity of 
the model. 



S 



The dependence of the self-energy functional Q on the fictitious deviation 5 is calculated 
at different chemical potentials fi which determines the average electron density of the lattice 
by the thermodynamical equation: in) = — |^ where N denotes the number of lattice sites. 
Fig. [2] is a plot of the difference Q(S) — Q(S = 0) as a function of S at U = 4 and T = 
for jj, = 2.0 (half-filling) and \x = 1.36 (off half-filling, n w 0.76). Notice that parameter 
5 can be varied only within the interval [0,1]. There are several stationary points in the 
two plots but, generally, the lowest point is usually the best guess for the proper stationary 
point representing physical properties of the ground state. The red arrows mark the proper 
stationary points in both cases respectively. At \i — 2.0 (half-filling), the proper stationary 
point is at 5 = implying that a uniform self-energy is preferred, while at fi = 1.36 (far 
below half-filling) , a non-zero fictitious deviation shows up which elucidates the possibility 
of a nematic state. Although it is not the exact stationary point of the whole variational 
space because other variational parameters are fixed, usually the true stationary point is 
not too far away, so the ground state at the true stationary point is very likely to remain 
nematic. 

We also use 2x3 clusters as the reference system to justify the results related to nematic- 
ity Fig. [3] also shows a plot of Q(5) — Q(5 = 0) as a function of 5 and nonzero stationary 
point which demonstrates the existence of a nematic ground state of the original lattice. 
This figure is very similar to Fig. [2] with stationary points at 5 = at half-filling and 5^0 
at n m 0.76, but the nonzero fictitious deviation 5 within this reference is smaller than that 
within the 2x2 reference system. As discussed in Sec. „ the parameter 5 measures the 
residual effect of nematic symmetry breaking from other parts of the whole lattice, so the 
larger the cluster, the smaller the residual effect, i.e., smaller 5. Therefore, 5 is expected to 
be smaller in the 2x3 cluster-based calculation. 

Antiferromagnetism and Locally Nematic State 

There are many competing phases in strongly correlated materials; whether two given 
phases can coexist is also an issue very much tied to the correlations of these phases. Usually, 
order parameters of two phases are calculated at the same electronic density n: if both order 
parameters are nonzero at the same n, the two phases coexist. Here, in order to determine 
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FIG. 2. Dependence of the grand potential Q on the fictitious deviation 5 is obtained from VCA 
calculation (the difference £1(5) — Q(5 = 0) is plotted) on an infinite square lattice tiled by 2 x 2 
clusters (see Fig. [I]). The on-site Coulomb interaction is U = 4 and temperature is T = 0. The 
energy unit is given by the nearest neighbor hopping of the original lattice. The calculation is 
looking for the stationary (minimum or maximum) point of £1(5) and the lowest stationary point is 
the best approximation to the ground state. The upper figure is calculated when chemical potential 
/i = 2.00 (half-filling) . The lowest non distorted stationary point with 5 = implies that there is 
no nematicity for half-filling ground state. In contrast, the lower plot at \x = 1.36 (off half-filling, 
n Rj 0.76) has a proper stationary point with 5^0 shown by the red arrow at 5 ~ 0.17. The 
nonzero 5 value suggests the nematic phase as a possible ground state solution in a single band 
Hubbard model. We further investigate 5 up to 1; in the region 0.36 < 5 < 1 both plots are 
monotonically decreasing so there are no other stationary points (not shown in this figure). 

the possibility of the coexistence of AF order and nematic order, the effective Weiss field M 
and the fictitious 5 are both processed in the variational calculation as discussed in Sec. „ 
We test the possibility of coexistence of the two states by comparing parameters M and 
5. Although they are not order parameters because they are not physical observables, they 
can still provide a "yes" or "no" answer to whether the specific state exists at a certain 
electronic density n. Therefore, as long as M and 5 are nonzero at the same n, the AF order 
and nematic order can coexist. 

After searching for possible stationary points in the variational space consisting of both 
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FIG. 3. 2 x 3 clusters (the inset figure) are used as reference system. The fictitious deviation 5 is 
obtained from VCA calculation. The difference S7(<5) — £1(5 = 0) is shown at ^ = 2.0 (n = 1.00) 
and jx = 1.36 (re ~ 0.76). At half-filling, the stationary point is still at 6 = indicating a normal 
state and a non-zero stationary point at 5 = 0.14 is found at /i = 1.36. This shows the stability 
of the local nematic ground state in the pseudogap region at sufficient doping n ~ 0.76 and the 
smaller 5 value is caused by the enlarged cluster size. 

nematic and AF phases, the fictitious deviation S and the effective AF Weiss field M are 
calculated simultaneously for different electron densities n and the results are plotted in 
Fig. HI The plot clearly shows that when one variational parameter is nonzero, the other 
one turns out to be zero. This result is also verified by choosing an AF ground state as 
the starting point of the variational calculation. Fig. 2] indicates that the proper physical 
state of this simple Hubbard model cannot support the coexistence of static AF order and 
nematicity, i.e., a static AF order with an AF vector Qaf — i^,^) cannot simultaneously 
break the local C4 symmetry. 

Experiments have shown that the nematicity can exist in the pseudogap region, where 
some local magnetism is still retained. First, it is possible that this magnetism is purely 
temperature driven. The reason we did not see such coexistence in the ground state is 
probably due to the fact that the static AF order with Qaf = (tt, tc) requires an electronic 
structure with a perfect C4 symmetry, while the nematic state breaks that symmetry. It 
is possible that the AF state coexists with local nematicity when it is not static or has a 
differently oriented AF vector. For example, a spin density wave phase or Qaf = ( 7r /2,0), 

11 



U=4, T=0, n=2.0, n=1.0 




joojoojooj 

OOiOOiOOi 

ooiooioo! 

000000 
000000 





i U=4 and T=0 

— ■— AF Weiss Field 
i — • — fictitious deviation 








• • • • • 



-0.2 2 



0.70 0.75 0.80 0.85 0.90 0.95 
n 



FIG. 4. The fictitious deviation 5 and AF Weiss field M are plotted simultaneously as a function 
of the electron density n at U = 4 and T = 0. M/0 means there is a static AF order with AF 
vector Qaf = (vr, 7r) and 5^0 indicates the nematic state appears, if they both are nonzero at 
the same n, that means the two states can coexist. In the plot, when 5 ^ 0, M is zero and when 
M 0, 5 is always zero. This plot indicates that the static AF phase with Qaf = (tTj tt) and the 
local nematic state cannot coexist. 



because this kind of AF phase does not strictly require the C4 symmetry, the two phases 
may have a chance to coexist. 



Strongly Correlated Case 



In order to study the case with strong electronic interactions, we consider the model at 
U = 8. Although the fictitious deviation S here is not a physical observable, it shows the 
response of the local electronic structure to nearby clusters within the nematic domain, so 
5, in some sense, represents the strength of the nematic order in the local nematic domain 
as long as it is calculated within the same cluster model. 

In Fig. [SJ values of S at different electronic densities n for U = 4 and U = 8 are shown. 
At U = 4, the nematic state can exist from the overdoped region to the optimally doped 
region. The value of 5 decreases with increasing electronic density, which implies that the 
nematicity is suppressed by local magnetic domains which begin to form as the electronic 
density becomes larger. At U = 8, 5 behaves the same way as that at U = 4 in the overdoped 
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FIG. 5. The fictitious deviation 5 from VCA calculation is plotted as a function of electronic 
density n for [7 = 4 and 8. For [7 = 4 (open square), 5 is decreasing with electronic density up to 
n < 0.8. For U = 8 (open circle), 5 behaves the same way as it does for [7 = 4 within the overdoped 
region but increases up to about one around optimal doping at U = 8, which is probably a signal 
of the transformation into the stripe phase. 



region. However, around optimal doping, it jumps to about 5 ~ 1 at U = 8. This abnormal 
behavior at U = 8 points to a picture where the local hopping in one direction becomes close 
to zero, which suggests that the two-dimensional lattice is reduced to quasi one-dimensional 
chains. This strong coupling phenomenon is very similar to the stripe phase with charge 
and spin modulations coupled with broken rotational and translational symmetries found 
in some strongly correlated high T c superconducting materials. Such a strong trend to form 
inhomogeneous pattern of holes and spins is often referred to as charge and spin phase 
separation instabilities. The same abnormality was also reported by Okamoto et q/.flO]. 
where they have given similar explanations. However, in our opinion, the stripe phase 
simultaneously breaks both rotational symmetry and translational symmetry, but the VCA 
model does not examine the breaking of translational symmetry. In addition, as mentioned 
in Sec. „ S is an artificial field, so it is important to verify the tendency to form the stripe 
phase by monitoring relevant order parameters. 
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Order Parameter of Local Nematicity 



In order to examine possible local nematicity, we look for a suitable order parameter. 
It is our belief that a cluster-based approach such as the VCA, where local correlations in 
the Hubbard model are treated exactly, is ideally suited for probing local order. Important 
insights can be obtained by examining a local order parameter which can be chosen as 
the difference between a physical observable along x and y directions such as the resistivity 
tensor jlj], spin correlation pj or some other structure factor. For example, STS experimental 
work by Lawler et al. 8[ measured the local density of states (LDOS) within each "sub-unit- 
cell" with a very high real-space resolution, which allowed them to successfully identify local 
nematicity. They defined a real-space order parameter to quantify real-space nematicity 
based on the LDOS at different atoms. 

Here, in searching for a locally nematic state in the two-dimensional Hubbard model, we 
define a nematic order parameter based on the zero-time, local density correlation function 
p which is defined as 

p(r, t = 0) = (n{r, t = + )h(r' = 0, t' = 0)) - (n{r, t = + )) 2 (5) 

where n(r,t) is the electronic density operator at position r and time t, i.e., h(r,t) = 
( c t ( r ! t + + ) c o-( r > t)) with summation of different spins cr. The equation can be simplified 
by Wick's theorem 

p(r, t = 0) = G(r, t = 0+)G(-r, t = (T) (6) 

where G(r, t) is the Green's function of the lattice at position r and time t. In order to 
calculate the Green's function and density correlation, Eq. [H] is rewritten using the cluster 
notations 

p(r, t = 0) = g J dudu'e^-~ k '^G(X, k, u)G{-X, ~k', u') (7) 



k,k' 



where L is the number of /c-points along a chosen direction in the calculation and G(X, k, oS) 
is the cluster Green's function with notations defined in the Fig. [TJ We plot the position 
dependence of the zero-time density correlation along the x direction in Fig. |6j The figure 
shows that as r goes out of the cluster scale, the density correlation p drops drastically, 
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FIG. 6. The zero-time density correlation p is calculated along the x direction and plotted as a 
function of x, i.e., p = p(x,t = 0). p has a peak at x = and it quickly drops to a relatively 
smaller value as x moves to 1 lattice site away from the origin. When x reaches 2 lattice sites away 
from the origin, which is also out of the cluster scale, p becomes almost 0. As x increases further, 
p remains close to 0. This plot shows that the zero-time density correlation is strongly localized in 
our model. The plots are calculated at U = 4, T = and electronic density n = 0.81 and a 10 x 10 
fixed mesh is used to sample the /c-points. 



which indicates that p is strongly localized within a small length scale. Except for the trivial 
correlation to itself, the most significant density correlation falls at the nearest neighbors. 
Let us denote the position of the nearest neighbor along the x direction as x and that along 
the y direction as y. Then, the nematic order O which breaks the C4 rotational symmetry 
can be defined by the difference between density correlations of the nearest neighboring sites 
x and y. 



p(x,t = 0)-p(y,t = 0) 
p(x,t = 0)+p(y,t = 0) 

In the normal (i.e., non-nematic) state, because of C4 symmetry, the local correlation func- 
tion, as defined above along the x and y is always zero. For a locally nematic state, the local 
correlation functions along the x and y directions will not be the same. Therefore, O is zero 
in the normal state and nonzero in the locally nematic state. Fig. [7] shows the evolution of 
the nematic order parameter O with electronic density n for U = 4 and 8. There is a clear 
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phase transition between a nematic phase and a simple static AF phase. Since there is no 
evidence for coexistence between these two phases in our model (as discussed in the Sec. 
the figure shows clear phase boundaries, which are at about n = 0.82 at U = 4 and n = 0.86 
at U = 8, between a C4 symmetry phase and a broken C4 symmetry phase correspondingly. 
The nematic order parameter O in the nematic phase does not change too much with elec- 
tronic density and indicates a relatively stable nematic phase. When compared with Fig. [5j 
where the fictitious deviation 5 jumps to nearly 1 at U — 8 within the region 0.8 < n < 0.86, 
O does not show similar abnormal behavior. Although 5 shows a tendency for the lattice to 
go into a stripe phase, the order parameter does not provide any support for the existence 
of this phase. In addition, the fictitious deviation field 5 is an artificial field and therefore, 
the possible existence of a stripe phase in the lattice is still an open question. 

The order parameter O is defined using the density correlation function p between nearest 
neighbor sites. As shown in Fig. The correlation function p has a peak at the origin 
and decreases drastically at distance 1, which suggests that the correlation function decays 
quickly as distance increases, so the electronic correlation is a short-range effect. Since 
electronic nematicity is caused by electronic correlation, the electronic nematicity is also 
a local effect. The value of p drops to zero beyond the scale of the cluster due to the 
fact that the VCA is a mean field treatment beyond the cluster dimension. The mean 
field approximation suppresses fluctuations and hence the correlation goes to zero. The 
nematic ground states are two-fold degenerate with different orientations. In the symmetric 
lattice, the degenerate nematic ground states spread locally over the lattice. The randomly 
disordered local nematic orders with different orientations cancel out each other, so the 
whole lattice, at large length scale, is symmetric. Therefore, the nematic fluctuations in the 
uniform Hamiltonian can exist locally within each "sub-unit-cell" and this conclusion gives 
strong support to the experimental results by Lawler et al. 

Dynamical Spin correlation function 

Magnetic neutron scattering is a very useful tool to investigate magnetic and density 
fluctuations of a given material. The magnetic neutron scattering spectrum is related to the 
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FIG. 7. Local nematic order parameter O is plotted as a function of electronic density n at U = 4 
(open squares) and U = 8 (open circles). The temperature for both plots is T = 0. This result is 
obtained by summing up to 10 x 10 momentum /c-points in the cluster reciprocal lattice. O is non- 
zero in the nematic phase with a broken C\ symmetry and is zero in the phase with a symmetric 
order. The system undergoes a phase transition at about n = 0.82 at U = 4 and n = 0.86 from a 
nematic phase to an AF phase. O in the nematic phase does not change too much as the electronic 
density n changes. 

dynamical spin correlation function. It is denned as 

S(r i ,r„t) = (S(r l ,t)S(r„0)) (8) 

1 1 r+oo 

S(fc^) = TfE^ / dte-^-^S(r t , rj ,t) (9) 
Jy —■ 2tc J-oo 

The spin correlation of the 2D lattice includes longitudinal S±(vi, rj, t) (S±(k, co)) and trans- 
verse components Snfri, rj, t) {S\\(k, u)). The two components of real-space correlation func- 
tions can be written as 

S±(r l ,r J ,t) = (S z (r t ,t)S z (r v 0)} 

= J2 G a(rj,ri,-t)G a (ri,r j: t) (10) 

a 

= l - (S + (r u t)S-(r j} 0) + S-fa, t)S + (r J; 0)) 

= lT,[G(r i ,r„t) + G(r J ,r l ,-t)} (11) 

Z a 
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where G a is the real-space-time Green's function and a is the spin index (up or down). The 
longitudinal correlation function is the same as the result of density correlation function. 
The momentum-space correlation functions can be obtained by Fourior transform of the 
real-space correlation functions. Because we perform most of the calculations within the 
cluster framework, the DSF is finally written using cluster parameters defined in Fig. [T] as 

S ± (k,u>) = f~ £ e"«*« £ G a (X jh q,u) 

J Z7FJV RBZ X ijt a 

xG a (Xij,k + q-Q,u + u) (12) 

S ll (k,L0) = W{ G ( k ^)+ G (- k i-")} 
^ a 

= ^ £ [e- lkXl] G a (X %]1 k, to) + e^G^Xij, -k, -co)} (13) 

Xij,a 

where Xij is a lattice vector between the cluster sites % and j; k, q and Q are reciprocal 
lattice vectors as defined in the caption of Fig. [TJ G is the Green's function in the cluster 
representation; N denotes the number of sites in the lattice. 

S±(k,oj) represented by Eq. IT2l involves coupling between two quasi-particles at q + k and 
q' with energy u + v and v' respectively in (k, u>) space. The correlation spectrum is not zero 
only when q + k = q' + Q and u + v = v' are both satisfied, i.e., v(k + q) — u(q) = u. This 
condition corresponds to an electron-hole excitation, which is only valid around the Fermi 
surface (k F ) of the cluster superlattice. S\\(k,u) is the momentum-space green's function of 
the lattice. 

The longitudinal and transverse spin correlation functions for magnetic excitations are 
calculated around in the momentum space at different energies. Fig. E] shows a slice cut 
of the DSF spectrum along (^,0) and (0,k y ). The large differences of the spectrum along 
x and y directions clearly show the asymmetry of the momentum-space density correlation 
function. This confirms that although the nematic state is a local effect in real-space, it still 
has a strong influence on some properties within the momentum space, especially on the 
two-particle correlations. 

Although our detailed theory only underlines the interplay of electron and spin density 
nematicities, the complete relationship between these order parameters and the high temper- 
ature superconductivity still needs further investigations. Strong competition of magnetism 
and superconductivity observed experimentally in both iron-pnictides and iron-chalcogenides 
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FIG. 8. The slice cut of the longitudinal and transverse spin correlation function S]_(k,cu) and 
Su(k,u) along (k x ,0) and (0, k y ) at U = 4 and T = with uj = 0.001 and electronic density 
n = 0.81. The result of S±(k,u) is obtained by summing up 20 x 20 momentum space points and 
integrating over the energy window from —1.0 to 1.0. Both plots show large differences along the 
x and y directions, which elucidates that the nematic properties. 

implies that a complete description of nematic ordering must involve also the superconduct- 
ing phase. The local electronic nematic ordering, which breaks C4 symmetry, can be coupled 
to the introduced nematic order parameter O by the fictitious parameter 5 through the next- 
nearest or next-next-nearest neighbor couplings. This is one of the directions that can be 
pursued for studies of the extended Hubbard model in the future. 

SUMMARY 

In conclusion, we have identified a local nematic phase away from half filling (in the 
overdoped region) by using the variational cluster approach to solve the one-band, two- 
dimensional Hubbard model. However, we do not find any coexistence of antiferromagnetic 
order and nematicity for their distinct symmetries. The variational cluster approximation 
shows evidence for electronic nematicity of the states consistent with the scanning tunnelling 
(spectroscopic-imaging) microscope measurements of the intra-unit-cell states in underdoped 
Bi 2 Sr 2 CaCu 2 0g + s [8J. The breaking of rotational symmetry by the electronic structure 
within each Cu0 2 unit cell becomes predominant in the pseudogap phase as the density of 
doped holes is reduced. Our results directly demonstrate that the nematicity can arise within 
a one-band isotropic (non-degenerate) Hubbard model. The nematic phase is suppressed as 
the electron density increases toward the half-filled state, but as the electron density gets 
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close to optimal doping, while the Coulomb interaction is large enough, the nematicity 
suddenly gets enhanced and is likely to give rise to a phase similar to the stripe structure. 
The dynamical structure factor shows that the real-space local effects also influence the 
correlations in the momentum space and the asymmetry of the correlation pattern changes 
with different incident energies. This report is one of the first studies where a locally 
nematic state in the single band Hubbard square lattice is obtained without any additional 
input (such as next-nearest neighbor couplings, inter-site Coulomb interaction or multi-band 
structure). 
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